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Abstract 

The shear-transformation-zone (STZ) theory of plastic deformation predicts 
that sufficiently soft, non-crystalline solids are linearly unstable against form- 
ing periodic arrays of microstructural shear bands. A limited nonlinear analy- 
sis indicates that this instability may be the mechanism responsible for strain 
softening in both constant-stress and constant-strain-rate experiments. The 
analysis presented here pertains only to one-dimensional banding patterns 
in two-dimensional systems, and only to very low temperatures. It uses the 
rudimentary form of the STZ theory in which there is only a single kind of 
zone rather than a distribution of them with a range of transformation rates. 
Nevertheless, the results are in qualitative agreement with essential features 
of the available experimental data. The nonlinear theory also implies that 
harder materials, which do not undergo a microstructural instability, may 
form isolated shear bands in weak regions or, perhaps, at points of concen- 
trated stress. 



I. INTRODUCTION 

In a recent publication M.L. Falk and I proposed a "shear-transformation- zone" 
(STZ) theory of plastic deformation in amorphous solids. This theory was based on direct 
observations of micro-scale rearrangements in molecular dynamics simulations. The central 
feature of the theory was the recognition that the STZ's — localized regions where irre- 
versible molecular rearrangements occur — are two-state systems; they can switch back and 
forth between only two orientations. Thus, the density of zones and the average of their 
orientations are appropriate order parameters for characterizing the internal states of such 
systems. The equations of motion for these order parameters, including terms that describe 
annihilation and creation of STZ's during deformation, produce a relatively simple theory of 
plasticity that exhibits strain hardening, a dynamic transition between viscoelasticity and 
viscoplasticity at a yield stress, and hysteretic effects. 

As presented in [^, the STZ theory is strictly a mean field theory; it allows no spatial 
correlations in the stress or strain fields or in the new order parameter. This assumption 
of spatial homogeneity, however, is inconsistent with the original picture when examined in 
closer detail. A localized shear transformation produces a quadrupolar elastic displacement 
field in its neighborhood, and a corresponding change in the local stress pattern. It is easy 
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to see that this change in stress may induce shear transformations in neighboring zones 
along the directions of maximum shear. (See Sections |T| and |I^ for details.) If each 
localized shear transformation induces others along certain preferred directions, then the 
system as a whole may be unstable against formation of a pattern of shear bands. This 
argument is the basis of numerical simulations described in an important series of papers by 
Bulatov and Argon. [0-^ The latter authors used a model of localized shear transformations 
that is more phenomenological than the STZ theory and that does not contain some of its 
central features, for example, the transition between viscoelasticity and viscoplasticity or 
the hysteretic effects. One purpose of the present investigation has been to explore the 
Bulatov-Argon picture of microstructural shear banding in the context of the STZ theory. 

In its original mean-field version, the STZ theory predicts the time dependence of the 
plastic shear strain e^'(t) in experiments (creep tests) where a uniform (deviatoric) stress 
s is applied quickly and then held at a constant value. For stresses less than an ultimate 
yield stress Sy, the strain approaches a value £j'„a;(s) = eP\t — >• oo) and remains constant 
there; that is, the two-state systems become "jammed." The quantity £^finaii^) diverges as 
s approaches Sy from below, as does the characteristic time for approach to the final state. 
A graph of s as a function of e^ji^ai therefore looks like an ordinary stress-strain curve with 
plastic yielding, although the STZ interpretation of it as the infinite-time limit is not strictly 
conventional. For applied stresses above Sy, the system comes to steady state at a nonzero 
strain rate s^jinaii^)' which increases with increasing s > Sy. In this case, the characteristic 
time for approach to steady state diverges as s approaches Sy from above. 

I have (so far) found only two examples of constant-stress experiments of this kind 
reported in the literature, an early one by Ender in 1968 and another more recent in- 
vestigation by Hasan and Boyce in 1995 0, both using the polymer PMMA. The work of 
Hasan and Boyce (HB) is rich in detail and is supplemented by a theoretical analysis of their 
data. At small applied shear stresses (up to about 80% of what I interpret to be an ultimate 
yield stress), the HB results seem to be consistent with the mean-field STZ prediction, i.e. 
the strain rate decreases monotonically to zero. At higher stresses, however, the curves of 
strain versus time have inflection points where the strain rate goes through a minimum and 
then increases, indicating the onset of some new mode of viscoplastic deformation. Signifi- 
cantly, the transition from viscoelastic to viscoplastic response occurs within a very narrow 
range of applied stresses; that is, the strain rate is highly sensitive to the stress within this 
range. I propose that this new mode of deformation may emerge from a microstructural 
shear-banding instability. 

The most common measurements of plastic deformation are performed under conditions 
of controlled strain, specifically, at constant strain rates. For soft materials such as polymeric 
glasses or soils, stress-strain curves generated in this way often exhibit strain softening. That 
is, the stress initially rises with increasing strain but then goes through a maximum and 
decreases, so that the stress ultimately required to sustain the fixed strain rate remains lower 
than the peak stress. Hasan and Boyce show clear examples of this kind of behavior; 
similar data for soils can be found in the book by Ishihara 0. The converse of the strong 
stress sensitivity seen in the constant-stress experiments is the fact that the peak stresses in 
the constant-strain-rate experiments are only weakly dependent on the strain rate. I suggest 
that this version of strain softening also can be understood as the result of an instability that 
sets in near the peak stress and produces a pattern of microstructural shear bands which 
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account for the deformation rate in the subsequent steady state. 

The theoretical analysis of Hasan and Boyce provides important background for the 
present investigation. Like the STZ theory, the HB theory is based on a picture of localized 
shear-transformation sites whose populations and states of deformation are the internal state 
variables that characterize the system. Unlike STZ, however, HB do not try to model the 
molecular rearrangements taking place at their sites or the interactions between these sites; 
thus they do not attempt to discover a connection between strain softening and shear local- 
ization. Nevertheless, there are deep relations between these theories. Both are attempts 
to identify the physically relevant internal variables whose equations of motion must ulti- 
mately be the basis of any satisfactory, history-dependent theory of plasticity — as opposed 
to, say, purely phenomenological theories based on yield surfaces p or Masing rules [§. 
(Unfortunately, HB use the strain itself as a state variable, which cannot be correct from a 
fundamental point of view.) Moreover, both theories identify cumulative internal changes 
during plastic deformation as the basis for dynamic yielding at high stress. 

The scheme of this paper is as follows. Section || contains a brief outline of the basic 
ideas of the STZ theory and serves to define concepts and notation for the rest of the anal- 
ysis. Some parts of the discussion presented there have not been published before now. In 
Section I derive an approximate, quasilinear, but fully tensorial generalization of the 
STZ equations that allows me to consider spatially inhomogeneous stress changes. Section 

contains a linear stability analysis of the spatially uniform steady-state solutions of the 
quasilinear theory. The expected instabilities emerge, but with some surprises. The quasi- 
linear analysis provides the basis for a nonlinear integro-differential equation that describes, 
in an effectively one- dimensional situation, how shear patterns may emerge in response to 
shear stresses. The solutions of the latter equation exhibit the strain-softening behaviors 
described in the preceding paragraphs. They also provide a mechanism for the formation 
of isolated shear bands in materials that do not exhibit strain softening. The nonlinear 
analysis, including numerical results, is described in Section |V|. 

The results presented here are, at best, a rough caricature of plastic deformation in shear- 
banding or strain-softening materials. Some features seem to be quite convincing, but others 
remain problematic. Answers to several central questions, especially the importance of the 
one-dimensional approximation made in Section 0, will require numerics and laboratory 
experiments that are beyond the scope of the present investigation. 

II. ELEMENTS OF THE STZ THEORY 

As in the discussion presented here pertains only to two-dimensional non-crystalline 
systems at effectively zero temperatures, i.e. at temperatures far below the glass transition. 
Also, following it is useful to start with the special case of a single preferred orientation 
for the zones and the applied stress, and to denote the number-density of zones oriented in 
the "-)-"/"_" directions by the symbol n±. Unlike however, the convention here is one in 
which the externally applied deviatoric stress is diagonal in the x, y coordinates, specifically 
= So, and s^y = 0. Choose the zones to be oriented (elongated) along the y 
axis, and the "— " zones along the x axis. With this convention, we can — temporarily — 
suppress tensor notation. 

As in Eq.(3.10), the plastic strain rate is: 
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= 6^ (Ae) W 



(2.1) 



where 6^ (previously denoted Vz) is the typical area of an STZ, (Ae) is the increment of 
local shear strain that occurs when one of these zones undergoes a transition (a number of 
order unity), and 

W = R+n+-R_n_. 



(2.2) 

" + " transitions respec- 
tively. These transition rates, at very low temperatures, are entropy-activated as opposed 
to thermally activated: 



Here, i?+ and i?_ are the rates for " + " 



- " and " - " 



R±{so) = Ro exp 



AV*{±so) 



(2.3) 



where AV^*(±so) is the excess free volume (area) needed to allow a "it" — > "ip" transition 
to occur, and f/ is the average free volume in the system. In |T[], Al^*(±so) was chosen to 
have the form 



Al^*(so) = V* exp(-so//i). 



(2.4) 



where V^* is of order the average volume per molecule and ft is generally much smaller than 
the shear modulus. Rq is an attempt frequency that, in [|l|, was found to be proportional 
to the square root of the total strain rate. Here, I shall assume that this factor is simply a 
constant. 

If there is no spatial dependence of any of the variables, then the equation of motion for 
n± is 



n± = R^n^ - R±n± + (e^' sq) A^. - Aann n± 



(2.5) 



This equation includes terms which describe creation and annihilation of STZ's at rates 
proportional to the rate at which plastic work is being done on the system, e^' sq, with 
coefficients Acr and Aann respectively. 

As previously, define tza = n_ — n+, ritot = n_ + n^. Also, 



R±{so) = Ro [C{so)±S{so)] 



(2.6) 



where 
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Then 



flA = 2RoS{so) ritot - 2RqC{sq) Ua - {e^^ So) Aann riA 



(2.7) 



(2.8) 



and 
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fltot = (^^' So) [2 Acr - Aann Utot] (2.9) 

Eq. (|2.9| ) implies that the ratio 2Acr/Aann = "^oo is the natural equilibrium value for 
ritot, the total density of STZ's. In principle, rioo could be a stress dependent quantity. 
Amorphous materials (even the Lennard- Jones solid in the molecular dynamics simulations 
described in |jl[]) generally exhibit a second order dilation in response to a shear stress; thus 
the equilibrium density of STZ's might increase as a function of Sq. For simplicity, however, 
I shall not consider that possibility here. 

Next, transform to dimensionless variables: n/\ = rioo iT-tot = "n-oo A, and write 

eP^ = ^w, (2.10) 

where 

M; = r(so)A-A; (2.11) 



and 



£o = ^ ( A5) ric 



T(so) = tanh 



^(so) 



2/?oC(so); 



^ smh £^ 



The equations of motion for A and A become 

w 



and 



A 



A 



1 -7SoA) 



7 



■wso{l-A). 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



(As) 



Here, 7 = (1/2)A„„„6- 

As pointed out in [0] (see also 
of solutions, one — the "jammed' 
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or 



00, which I shall assume to be a stress-independent quantity. 

11|] and H), Eqs. fpl^) and (|2T5|) have two kinds 
viscoelastic" solution — in which the plastic rate 
of deformation e^' vanishes in steady state, and another — the "viscoplastic" solution — 
in which it does not. In the first case, there is a one-parameter family of solutions with 
A = Aq = T(so) Ao where A = Aq can have a continuous range of values determined 
dynamically by the density of STZ's in the initial state of the system. In the second, 
viscoplastic case, we can have only A = Aq = 1/(7^0) and A = Aq = 1. The jammed 
states are dynamically stable attractors for values of sq less than an ultimate yield stress Sy, 
determined by 



1. 



(2.16) 



Conversely, the viscoplastic state is stable for sq > Sy. To see these stability properties, 
write A = Aq + Ai(t), A = Ao + Ai(t), and compute the amplification rate uJo{so) for which 
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Ai = uqAi, Ai = cuqAi. The result is that, for sq < Sy, 
(because Aq is undetermined), and ujq = Qq, where 



the two eigenvalues are 0^0 = 



fio(so) = 2C(so) 7SoT(so)-l 



(2.17) 



For So > Sy, the Ai and Ai modes are decoupled, and both have the eigenvalue ti^o = —flo- 
Because flo is an increasing function of sq that vanishes at Sy, the stable attractor crosses 
from the viscoelastic to the viscoplastic branch of steady states at sq = s. 



III. QUASILINEAR THEORY OF INHOMOGENEOUS STRESS PATTERNS 

So far, this analysis has pertained only to systems in which the shear stress has magnitude 
So everywhere and has only a fixed orientation. To incorporate inhomogeneous stress tensors 
Sij with varying orientations into this analysis, we need — at least temporarily — a fully 
tensorial version of the theory. No such version of the complete STZ theory exists at present. 
(See Appendix A of Falk's dissertation for a first attempt at developing such a version.) 
Therefore, I shall resort to a quasilinear approximation p!0| , p!2[] in order to derive equations 
of motion for the inhomogeneous stress variations. 

Start by noting that the (deviatoric) stress and plastic strain must in general be traceless, 
symmetric tensors, Sij and efj respectively. It follows from the form of Eqs.( p.lC ), ( p. 14 ), 
and (|2.15| ) that A, w, and T also must generalize to traceless, symmetric tensors, and that 
A must be a scalar. One way to obtain such a structure is to write 



T ^ 



Xsij, \^Vo/{vfofi); 



(3.1) 



that is, linearize the right-hand side of Eq.( p.l3D in sq and then let sq become Sij. To be 
consistent, we also must approximate the time constant t(so) by its limit at small so, which 
I denote simply by To: 



To 



(3.2) 



These estimates for A and tq, however, should not be taken literally; they represent only 
one of possibly many ways in which we might start from atomistic considerations to derive 
a theory with the general phenomenological structure shown below. 
With the approximations (|3.1| ) and (|3.2|) , we now have: 



■pi _ £o 

~ To 



W 



Wi 



A Sij A — Ajj. 



(3.3) 



Also, (using summation convention), 

1 



1"0 



Wij - -liwkiSki) Aij 



(3.4) 



and 



A 



2 To 



(wki sh)7(1 - A). 



(3.5) 
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The dissipation rate, denoted by sq in Section |I|, generalizes here to (1/2) e^' Sfc^. Note 
that, if we guess that A ^ 1 quickly and stably and therefore set A = 1 in Eg. (|3.3D, we obtain 



the quasilinear STZ equations that were used in and |ri|. If we did this here, however, 
we would miss some of the most important features of the shear-banding mechanism. 

The new element of the present theory is that the stress tensor Sij is to become a 
dynamical quantity varying with both position and time (x, y, and t). Before writing 
equations of motion for Sjj, I must specify the physical picture that I have in mind. The 
molecular dynamics simulations reported in imply that the STZ's are sparsely distributed 
sites at which small clusters of atoms undergo irreversible rearrangements in response to 
applied stresses. To visualize an STZ, think of a void in an elastic material, and place 
a small group of atoms inside this void in such a way that their average free volume is 
somewhat larger than that for most other atoms in the system. The void has some degree of 
structural stability; it can deform elastically but, because of the configuration of atoms on 
its surface, it resists collapse. Rearrangements of the atoms that are caged within the void 
couple to its shape and, therefore, to the stress field in the elastic medium in which the void 
is embedded. (The stress field outside the void is similar to that in the neighborhood of an 
Eshelby inclusion. |T^) This stress field — outside the STZ's — provides the mechanism by 



which the STZ's are coupled to each other. 

To distinguish it from the uniform applied stress Sq (as defined at the beginning of Section 
fl|), denote this fluctuating part of the stress field by the symbol Sij{x,y,t). For simphcity, 
first consider a "— " "+" STZ transition located at the origin of the coordinate system and 
aligned along the original x, y axes; and compute the stresses and strains in its not-too-close 
neighborhood by the potential method. That is, define ^{x,y) such that 



V2$ = 0; = 9,$; Uy = dy^; (3.6) 



and choose 



^, . y'^ - x"^ cos26 , , 

[x^ + y^Y 

so that Ux ~ x~^ along the x axis and Uy ~ —y^^ along the y axis. Then the deviatoric 
stress and strain fields have the form 



Sxx - Syy ~ 9,$ - dy^ (x^+y2Y = ^3.8) 



where r and 6 are the conventional polar coordinates. This function has negative lobes along 
the X and y axes and positive lobes along the lines x = ±y, so that the shear in the direction 
of the initial deformation is largest along the diagonals. 

Now suppose that the STZ transformation which caused the stress pattern ( |3.8| ) is rotated 
by an angle relative to the x, y coordinate system. Then 

(l>{x,y) = ^[{y^ -x^)cos2^-2xy sm2^] = -^Vij{^), (3.9) 
where V is proportional to the director matrix: 

V{^) = constant x f ""^^^t ^'""^fo) • (3-10) 
\sm2v — cos 2v J 
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If, according to ( p.6|) , $ is the potential for displacements, then, for any position r far 
enough away from the origin that it is outside the transforming zone, 

4(r) = 2^e^iv) = -2/z6^a.9, ^4(0). (3.11) 

Here, efj{r) is the elastic strain rate induced at r, /i is the elastic shear modulus, and b is 
the characteristic size of an STZ. The quantity i^ii^) strain rate induced by the shear 

transformation at distances r of order b, which, in turn, is proportional to the rate nA(0) at 
which the STZ population is changing at the point r = 0. More generally, the source of the 
fluctuating stress field Sjj(r) has a density distribution at a point, say, r' given by 

ntotSMir') = (Ae)nA(r') = (As) Akiir'). (3.12) 

Integrating over the population of STZ's with this source density, we find an equation of 
motion for the deviatoric stress field: 



Sij{v) = -2^b'^ {Ae) Uoo j (f'r' 



d d (vk - r'f,) {ri - r'l) 
dvi dvj |r — r'l^ 



Akiir')Hi\r-r'\-b). (3.13) 



Here, H is the Heaviside function that cuts off the integration inside the core of the STZ. 
(It may be more realistic, but probably not essential, to use a smoother cutoff.) Because 
this stress field comes directly from Eqs.( |3.6| ) and ( p.7|) , it automatically satisfies the force 
balance and compatibility conditions. 

Eq.( |3.13| ) is best expressed in terms of the Fourier transform of %(x, y), say Sjj(k). Then, 

2,iAlwn^ = ^^('^^ + ^^^"^^ 

ky{k) = A,M - ^i(k) A,,(k). (3.14) 



2fi (Ae) 62 n. 

Note that these are linear relations between the time derivatives of the fluctuating part of 
the stress field and the order parameter A. 

The kernels Q contain much of the dynamical information that we need. Write the 
wavevector k in the form = kcosip, ky = A; sin?/'. Then 

Giik, ^/j) = -cos{A^) g{bk); g2ik, ^) = - sm{A^) g{bk); (3.15) 

and 

roo dp r2TT 



Jk p-^ Jo 



Jo 

= 127r/€M -^J4p= Mk), 3.16 

Jk p-^ K 

where k = bk and the J's are Bessel functions. Therefore, g{K,) approaches ttk'^/A as k — > 
and is proportional to k"^^'^ as k — oo. It rises through a maximum of about 4.176 at 
K = Kc — 3.611 and then oscillates as it decreases. In the following, I denote the maximum 
value of ^(k) by z/~^, i.e. z/c = 0.2395. 
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IV. STABILITY OF THE UNIFORM STEADY STATES IN THE QUASILINEAR 

THEORY 



The next step is use the results of the preceding Section to examine the hnear stabihty 
of spatially uniform steady-state solutions of our equations of motion. As we shall see in 
Section 0, the linear stability analysis of steady states misses crucially important transient 
phenomena, but it does give us some insight about the behavior of this system. 

It is convenient to use a notation in which 

% =fn° ^ ) + e^"-'""'^*, (4-1) 



with analogous notation for components of A, and 



A = Ao + Ai e*'^-'-+^*. (4.2) 



The second term on the right-hand side of Eq. ([4.1|) is the fluctuating stress tensor Sij. Like 
all other first order terms in this linear analysis, its time dependence is exp(co't). 
At zero'th order, according to Eqs. (|3.3|), ( |3.4| ), and ( |3.5| ), we have 

Ao = -(AsoAo-Ao)(l-7SoAo), (4.3) 

To 

and 

Ao = ^ (Aso Ao - Ao) (1 - Ao). (4.4) 

To 

In this quasilinear approximation, the ultimate yield stress Sy, determined by Eq. (|2.16| ), is 

Sy = I/a/tA- 

Consider first the jammed situation for which sq < Sy, and examine the linear stability 
of the stationary state with Ao = AsoAo. Aq remains a free parameter to be determined by 
initial conditions. (We shall compute Aq explicitly in Section 0.) 

The first-order equations are the following: 

w To Ai = (A Ao si + A So Ai - Ai) (1 - 7 sq Ao); (4.5) 



wtqAs = AA0S2 - A2; (4.6) 

and 

w To Ai = 7 So (A Ao si + A so Ai - Ai) (1 - Ao). (4.7) 

From these relations, we find that 

XAoDa{so) . " AAo . , . 

Ai = , ^ ^1' ^2 = — S2; 4.8 

where 
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(4.9) 

The combination of Eqs.( |3.1^ ) and ([4.8|) is a pair of hnear, homogeneous equations that 
determines the amphfication rate uj{k, ip): 

s^ = -cos{4^)u^g{bk) ( ^^^'nf\ l - sin(4V^) z/^o ^(&A:) (^^) S2; (4.10) 

\ujTo + D{so)l Vu;ro + l/ 



and 



S2 = - sin(4^) m) ( ^^^"nf\ l + ^^^(4^) ^(&A;) f ^2. (4.11) 



Here, 



u^ = 2fx\{AE)b^n^. (4.12) 



As we shall see, u^o is the dimensionless group of parameters that controls stability against 
microstructural shear banding. 

To understand the nature of the stability spectrum, it is easiest to look at two special 
cases. The most important situation for our purposes is the orientation ip = ±7r/4, where 
the bands are aligned with the direction of maximum applied shear stress. In this case, the 
si and §2 modes decouple, and we have: 

u{k,n/A)ro = ( ^(^^) Ao^a(^o) - /^(^o) {h mode) 3 
[-Uoog{bk)Ao-l (s2 mode). 

As might be expected, an instability (positive u) can occur in the §1 mode where the shear 
deformation is parallel to the bands; the S2 mode, where the shear is normal to the bands, 
always decays. Suppose, for the moment, that Aq ^ 1 during the initial transient — a 
possible consequence of the dynamics described by Eq.( [4.4| ). Then, for the Si mode in 
Eq.(ra, 



euro 



si 



v^g{hk)-l. (4.14) 



There is an instability for all values of Sq < Sy, but only for u^o > I'c — 0.2395. When 
the latter condition is satisfied, the instability occurs in some range of values of bk in 
the neighborhood of Kc — 3.611. (The numerical constants Uc and Kc are defined following 
Eq.( pj.l6D .) Note that the first factor on the right-hand side of Eq.( |4.14] ) vanishes as Sq ^ Sy, 
implying that the instability weakens at large stresses and disappears at the yield stress. 

Unexpectedly (for me), the instability is stronger in the orientation ip = 0, where the 
si and §2 modes again decouple, but the bands are at forty-five degrees from the direction 
of maximum applied shear stress, i.e. in the direction in which the applied shear stress 
vanishes. We have 

^ ^ \uoogibk)Ao-l (s2 mode). ^ ' 
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In this case, only the §2 mode can be unstable. The shear deformation is again parallel to 
the bands, but there is no weakening factor analogous to that in Eq. (|4.14 ) for this mode, 
and therefore it will grow more rapidly. Apparently, at ip = n/4, a. substantial fraction 
of the STZ's are already aligned with the applied stress by the time that the system has 
reached its uniform steady state, and thus fewer zones are available to be transformed by 
the induced stress si. No such suppression occurs ai ip = 0, because there the bias in the 
orientation of the zones is not aligned with the orientation of the emerging bands. 

To complete this analysis, look at the viscoplastic states with sq > Sy. In steady state 
with nonzero plastic strain rate, Aq = 1/(7,50) and, without ambiguity, Aq = 1. In place of 
Eqs.(|4.8|), we have: 



Ai 



WTo + D{sq) 



Si, 



A, 



A 



S2, 



where 



(4.16) 



^(.0) = § - 1. 



(4.17) 



The analogs of Eqs. (|4.1C ) and ( 4.11|) are 



,x N / D{so) (Sy SoY\ . 
si = cos 4V^ 1^00 g{bk) , . si 



sin(4V^) Uoog{bk) 



1 



UTo + l + D{so) 



S2] 



and 



S2 



V UJTo + D[Sq) J 



S2- 



From these we find, for = 0, 



cc;(A;,0)ro 



D{so) \v^g{hk) {sy/s^f - 
1 + D{sq)\ [z/oo g{hk) {sy/sof - 1 



and, for ip = 7r/4, 



u;(fc,7r/4) To 



Diso) [u^g{bk) {sy/soY + 1 
l + D{so) yoog{hk){sy/so 



+ 1 



(si mode) 
(s2 mode); 

(si mode) 
(s2 mode); 



(4.18) 



(4.19) 



(4.20) 



(4.21) 



Now both modes are unstable at ^/^ = 0, and both are stable at ip = 7r/4. The instability 
occurs only in a finite range of stresses 1 < [so/syY < {i^oo/^c)', thus, it can be suppressed 
by driving the system further into the viscoplastic regime by increasing the applied stress 
Sq. Note that the geometry of this situation is similar to that of the crazing instability in 
polymers. The pattern of most rapidly growing modes is aligned with the principal axes of 
the stress tensor rather than along the direction of maximum shear. 
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V. NONLINEAR ANALYSIS 



We must look at nonlinear behavior in order to understand the implications of these 
linear instabilities. To make progress along these lines without a great deal of new formal 
development, we can restrict ourselves to situations in which the principal axes of the stress 
tensor retain a single, fixed orientation everywhere throughout the system and at all times. 
Within this restriction, we can examine the Si mode that is unstable for sq < Sy, because the 
applied and induced stresses in that case remain everywhere aligned along the x, y axes, and 
all spatial variations occur only along the diagonal orientations ip = ±7r/4. While this mode 
is not the most strongly unstable, it is the one that couples to the applied stress and thus 
describes the deformations that would be measured in response to the applied loads. The 
restriction of fixed stress orientation substantially limits the conclusions that we can draw 
from this analysis. On the other hand, there is no longer any reason to use the quasilinear 
approximation. Several physically important features of the following results emerge from 
the strong stress dependence of the full STZ theory. 

It is convenient to make the following changes of variables 



"yy 



7 



-A 



yy 



A(e,t), 



(5.1) 



and ^ = (x ± y)/b\/2; i.e. ^ is the spatial variable in units b along one or the other of the 
diagonal directions. For simplicity, express time in units Rq^. With these transformations 
and the additional spatial dependence, Eqs. (|2.14|) and (|2.15|) become 



A = 2C(a) [T{a)A-A] (1 - a A) 



(5.2) 



and 



A = 2 a C{a) [T(a) A - A] (1 - A). 



(5.3) 



Here, 



C(cr) = exp 



—a cosh 



a 



cosh 



a sinh 



a 



(5.4) 



tanh 



a sinh 



-) 

0-1/ 



(5.5) 



and a = Vq/vj, a\ = 7/i. Symbols such as C((j) now denote functionals C[(j(^,t)]. The 
yield stress 0"^, according to Eq. (|2.16|) , is determined by 



^y'^^^y) = 1; 



(5.6) 



Tightly packed, hard-to-deform materials such as glasses at low temperatures or the 
Lennard- Jones material studied in |l[] have small free volumes Vf and therefore large values 
of a. Conversely, softer materials such as polymers or clays should have smaller values of 
a. For large a, T{a) = 1 at all values of a appreciably larger than ai/a. Thus, Eq.( p.6|) 
imphes that ay is of order unity or, equivalently, Sy is of order 7"^. In [||, the value of Jx 
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was found to be of order the yield stress, roughly two orders of magnitude smaller than the 
elastic shear modulus. This estimate, along with the preceding estimate for 7, then implies 
that (Ji = 1. 

To be consistent with previous notation, write cr(^,t) = cro(t) + 5"(^,t), where ao(t) is 
the applied stress and a is the induced stress fluctuation. Because only the si mode plays 
a role in this situation, we can write Eq.( p.l3| ) in the form 



a = v g * A, (5.7) 

where u = 2/17 (Ae) 6^ rioo- The quantity 2/^7 is approximately the ratio of the shear 
modulus to the yield stress, and is therefore very roughly of order 10^. The quantity (Ae) 
is of order unity. As seen in the stability analysis, we need values of u greater than about 
0.25 to generate instabilities. Accordingly, 6^ rioo must be of order 10~^, which is consistent 
with our basic assumption that the STZ's are widely separated from each other. 

In Eq. ( |5.7|) , the symbol g* denotes the integral operator obtained by Fourier transforming 

g * A(0 = Jd^' m e'^^^-^'^ Ha- (5.8) 
Using Eq.( |3.16|) , we find 

g*m = j\'g{C-OHO (5.9) 



where 



g{^) = — L [3 sin(2^i) + 3 sin(4^i) + sin(6^i)] ; 9i = arccos^. (5.10) 



2e 

The function g{^) is shown in Fig. 1. It is strictly local; that is, it vanishes for |^| > 1 (an 
artifact of the sharp cutoff in Eq. (|3.13| )). The combination of Eqs. (|5l^ ), (|5.3| ), and ( ^.7|) is a 
nonlinear system of equations that can be solved for A(^,t), A(^,t), and a{^,t). 

Before starting this nonlinear analysis, however, note that the stability spectrum for the 
uniform state of the full STZ model is qualitatively different from the quasilinear result. 



Eq.( 4.13 ), in the case of hard materials with large a. For applied stress o"o < o'y, and for 
modes of dimensionless wavenumber k oriented along one of the diagonals, the amplification 
rate is 



2C((To) \l-aoAoT{ao)]Aor{ao)ug{K)-\l-aoT{ao)] . (5.11) 



This formula has the same structure as the Si mode in Eq. ( [4.131) . However, the quantity 
v g{K) is multiplied here by an extra factor T'(cro), which becomes exponentially small for 
c"o > cTi/a. Therefore, for large a and for ao not too small, the system is stable except for 
unphysically large values of D. The "weakening" factors in Eq. (|5.11|) (in square brackets) 



also slow the instability at stresses near cxy. At small values of ctq, on the other hand, the 
factor 2C((To) becomes small of order exp(— a), causing any instability that might occur 
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there to grow very slowly. We therefore expect to find little or no microstructural shear- 
banding instability in hard materials. We shall see, however, that the absence of instability 
does not imply the absence of shear banding. 

Turn now to the numerical solutions of Eqs. (|5.2| ), ( pl3|) , and ( jSTTl) . To solve these equa- 
tions, I have used a simple discretization of the integral kernel g{^), and have applied periodic 
boundary conditions, usually choosing the length of the system to be an integral number of 
wavelengths 2tt / k.^. of the most rapidly growing mode. In this way, the problem is reduced 
to a large set of coupled, ordinary, nonlinear differential equations. The initial conditions 
are always cr(^, 0) = and A(^, 0) = Aj„. We shall see that Aj„ provides crucial information 
about the way in which the system was prepared, specifically, the degree to which it was 
annealed. In order to allow an instability to develop, I have chosen A(^, 0) to be some 
small-amplitude, random white noise, corresponding to uncorrelated inhomogeneities in the 
initial orientations of the STZ's. Specifically, 



where the angular brackets denote a statistical average and M is the noise strength. 

For reasons mentioned previously, which I shall discuss further in Section |VJ it is not 
sensible to try to make quantitative comparisons between the results of this highly simplified 
calculation and real experimental data. Nevertheless, it is useful to use the results of Hasan 
and Boyce (HB) as a qualitative guide. First consider their creep tests in which a uniform 
applied stress cxo is turned on suddenly at t = and held constant. I assume that HB were 
using a relatively soft material, and therefore choose a = ai = 1. For these parameters, 
the yield stress is ay = 1.132. In their Fig. 4, HB show creep-test data for a sample that 
underwent a sharp transition from mean-field-STZ "jammed" behavior to rapid yielding 
over a narrow range of applied stresses, approximately 90 - 100 MPa. From this behavior, 
I guess that their ultimate yield stress — their analog of ay — was about 100 MPa, and 
that I should look at dimensionless stresses ao in some range near unity in order to make 
comparisons with their results. 

The sample used in this particular HB experiment was annealed so that the initial pop- 
ulation of transformation zones must have been quite small. It follows that we should look 
at small values of Aj„; that is, as in HB, we expect that the density of active sites is small at 
first and then increases as the system undergoes plastic deformation. Our stability analysis 
already gives us useful information about how this may happen. In Fig. 2, I have used 
Eq.( [5.1lD to compute the amplification rate for the most rapidly growing wavelength, u{kc), 
as a function of the stress ao- For lack of a better estimate, I have used u = 0.4 here 
and throughout the rest of this presentation, because that value produces moderate insta- 
bilities and, as mentioned previously, is consistent with rough estimates of its constituent 
parameters. To compute Aq, I have used Eqs.( pl^) and (|5.3|) to find 



The two functions uj{nc, ctq) shown in the figure correspond respectively to a highly annealed 
specimen with A^ = 0.1, and to one that is very rapidly quenched so that Aj„ = 1.0. The 
quenched specimen is strongly unstable at small stresses and stabihzes only just below ay. In 
contrast, the annealed specimen is unstable only in a narrow range of stresses from ao = 0.97 




(5.12) 



Ao = lim A(t) 




(5.13) 
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up to the yield stress, cxy = 1.132. At the smaller stresses, the density of zones never grows 
large enough to trigger the instability. 

The corresponding results of the nonlinear calculation for an annealed system with Aj„ = 
0.1 are shown in Fig. 3. Here, the plastic strain e^\t), in units Eq = 6^ (Ae) noo/2, is shown 
as a function of time (in units Rq^) for a sequence of applied stresses corresponding roughly 
to the stress range shown in HB Fig. 4. Specifically, the stresses shown in Fig. 3 are 
do = 0.4, 0.8, 0.9, 0.95, 0.975, 1.00, 1.05, 1.1, and 1.2. Again, i) = 0.4. The noise strength is 
M = 10~^. As expected, for stresses (Tq < 0.95, the strain rises to its uniform, steady-state 
value at which the system remains "jammed" indefinitely. There is a rapid transition in 
the range 0.95 < ctq < Cy, where the strain rate first slows but then increases to a new 
steady-state value, indicating the onset of microstructural shear banding. In this range of 
stresses, the system enters a stable flowing regime in which the dimensionless shear rate 



is localized in bands, as shown in Fig. 4 for the case ctq = 1.0 and time t = 100. As the 
stresses (Tq rise above ay, the banding instability becomes transient and the system reverts 
at long times to its mean-field STZ behavior for the viscoplastic regime. 

Now consider the second kind of experiment mentioned in the Introduction, i.e. the 
conventional situation in which a stress-strain curve is measured at constant strain rate. 
Assume that the spatially averaged total strain rate, say Eq Qq, is held fixed. Assume further 
that the total strain is a simple sum of elastic and plastic parts, and that the elastic strain 
e'^^ is given by Hooke's law, i.e. cr = 2 fi'j e^K Because the spatial average of the fluctuating 
stress a vanishes, the condition of constant go becomes an equation of motion for the applied 
stress: 



where q{t) denotes the spatial average of the plastic strain rate q{C,,t) defined in Eq. (|5.14| ). 

Figure 5 shows stress-strain curves computed using Eqs.(p.2|), (p.3|), and ( |5.7| ), supple- 
mented by Eq.( ^.15| ). The material parameters are the same as in Fig. 3, in particular, 
Ain = 0.1. The strain rates are: go = 0.008, 0.016, 0.032, 0.048, and 0.064. The peak 
stresses seen here are roughly comparable to those shown for the analogous HB experi- 
ments, and are in the transition range seen in the constant-stress calculations with a slight 
overshoot because the apphed stress continues to rise during the onset of the instability. 
The subsequent stress drops, however, are substantially larger and sharper than in HB. I 
have not so far found any way of correcting this discrepancy simply by adjusting parameters 
in the present theory, and therefore suspect that the problem lies with one or more of the 
major simplifying assumptions that are built into it. 

The behaviors shown in Figs. 3 and 5 represent only a very small part of the space 
spanned by the parameters a, u, Aj„, etc.. For example, the stability spectrum shown in 
Fig. 2 implies that creep tests performed on rapidly quenched soft materials, with large 
Ain, would exhibit shear banding instabilities at arbitrarily small applied stresses, and that 
these instabilities would actually weaken and disappear at larger stresses. Similarly, such 
materials (if they exist) would exhibit strain softening at constant strain rate only at low 
rates and with small peak stresses. It will be interesting to learn whether such behaviors, or 
others that may be less apparent results of this nonlinear model, actually occur in nature. 




(5.14) 



(5.15) 
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The last question that I shall address with these techniques is whether there may be some 
connection between the microstructural shear banding discussed so far and the formation of 
macroscopic shear bands in materials that do not soften in this manner. Consider a case in 
which the system is strongly stable against formation of microstructural shear bands, but in 
which there is a weak spot characterized by a locally high density of STZ's. Further, let the 
orientations of the zones in this weak spot be such that they will transform easily in response 
to an applied stress, and thus will induce additional stresses in their neighborhoods. 

The behavior of such a system is shown in Fig. 6. Here, as before, cxi = 1, and z> = 0.4, 
but a = 10. Initially, A(^,0) = and A(^,0) = 0.5 everywhere except in a narrow region 
(9 of 420 discretization intervals) at ,^ = 3.5 where A(,^,0) = —1, A(,^,0) = 2. There is no 
other heterogeneity and no noise in this initial state. A constant, uniform stress ctq = 0.5 is 
applied abruptly at t = 0. In the absence of the weak spot, this system would relax stably 
to a uniform, jammed state with g = after about 2000 time units. (The relaxation is 
slower than in previous examples because of the large value of a and therefore small value 
of C((To).) Figure 6 shows the resulting dimensionless strain rate t) at times t = 200 and 
t = 5000. At the earlier time, a central shear band and two weaker side bands are emerging, 
and the whole system is also deforming as indicated by the nonzero values of q well away 
from the bands. At the later time, all of the shear is taking place in the bands, and the 
peak strain rate in the central band is now an order of magnitude larger than it was earlier. 
This configuration seems to be completely stable; I have continued the calculation out to 
t = 10, 000 and find no further changes. 

This numerical experiment is as close as I can come to simulating a shear band in this one- 
dimensional truncation of the plasticity problem. In a real situation, the band would start 
at a stress concentration near some defect or surface irregularity, and would propagate away 
from that point along a direction of large shear stress. (See, for example, observations by 



Kramer JTSl.) Such intrinsically two or three dimensional behavior has no analog here. The 



most that can be said so far is that this calculation indicates the possibility of isolated band 
formation even in hard materials where no extended, microstructural shear-band patterns 
can occur. It also suggests that a macroscopic shear band might consist of a cluster of 
microstructural bands. 



VI. CONCLUDING REMARKS 

The STZ theory, including its extension discussed here, exhibits a wide range of be- 
haviors which look qualitatively like phenomena that occur in real materials. In contrast 
to other more phenomenological theories of plasticity, it has the advantage of being based 
on a microscopic picture of the internal states of noncrystalline solids; thus, it ultimately 
should be possible to compute the STZ constitutive parameters from first principles. More 
important, in my opinion, is the fact that the STZ theory consists of a small set of equations 
of motion for internal state variables, and therefore is intrinsically simpler and more general 
than conventional theories PJT6[]. The structure of the STZ theory, by itself and without 



reference to specific mechanisms, has experimental implications. For example, the diverging 
relaxation time near the ultimate yield stress in non-softening materials is a very general 
prediction that ought to be checked experimentally. Another possibility, as I have argued in 
a recent publication, is that the new dynamical degrees of freedom in the STZ model may 
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resolve the puzzle of how breaking stresses can be transmitted through plastic zones at the 
tips of brittle cracks. [|l^ 

Microstructural shear banding, while not so general a prediction as the diverging re- 
laxation time near the yield stress, seems to be an inevitable consequence of the specific 
STZ mechanism; the banding instability comes directly from the quadrupolar sjTiimetry of 
the STZ transformation. However, the theoretical analysis presented here is still far from 
definitive. Its most serious shortcoming is that it is only one dimensional in the places where 
it tries to relate to experimental observations, i.e. in Section 0. A more likely picture is 
the one shown in Fig. 6 of Bulatov and Argon where the two-dimensional pattern of 
shear bands consists of finite-sized patches oriented along the two equivalent directions of 
maximum shear stress. The way in which these patches interact and interfere with each 
other may be a crucial aspect of the pattern-forming mechanism, in which fully two- 

dimensional analysis will be necessary in order to understand the shear-banding instability 
in the STZ theory. 

There are several other serious shortcomings of the present analysis, most — or all — 
of which I believe can be remedied by plausible extensions of the existing theory. One such 
shortcoming is that the STZ theory assumes that there is only a single kind of zone rather 
than a distribution of them with a range of transformation rates. Hasan and Boyce 
include such a distribution in their constitutive model, and they use its width as one of their 
adjustable parameters for fitting experimental data. I suspect that the implicit assumption 
of an infinitely sharp distribution in the present analysis may be the principal reason why 
the stress drops in Fig. 5 are so abrupt and deep. It would not be difficult to rewrite the 
present theory as an average over different varieties of STZ's and thus smooth out such 
unrealistically sharp features. However, it would be much more satisfactory to derive such 
distributions from a first-principles theory of deformations in amorphous materials rather 
than just to introduce them phenomenologically. 

Another shortcoming is that I have not, so far, extended the STZ theory to finite tem- 
peratures. Including thermal activation in the rate factors seems a simple matter, and the 
temperature dependence of the phenomena studied here might provide useful tests of the 
basic theoretical ideas. 

Finally, there is the question of length scales in this theory. It seems puzzling that the 
only length scale to emerge in this analysis is b, the size of the zones. Why should the spacing 
of the microstructural shear bands scale simply with b rather than depending also on the 
density n^t? Why, for that matter, should such an instability occur only on microscopic 
rather than macroscopic length scales? Perhaps a better theory will produce a screening 
length for the interactions between zones. And perhaps, in a two-dimensional version of the 
theory, interactions between patches of different orientations will produce a characteristic 
patch size. There seems to be a wide range of possibilities both for extending the STZ 
picture and for testing its general validity. 
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FIG. 2. Amplification rate lo{kc) defined in Eq.( 5.1l| ), as a function of applied stress do, for two 
different values of according to Eq.( 5.13| ). The solid curve, for Aj„ = 0.1, corresponds to an 
annealed system with few defects; the dashed curve, for Ai„ = 1.0, corresponds to a system that 
has been quenched so rapidly that the initial STZ density is high. 
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FIG. 3. Plastic strain as a function of time for a sequence of simulated constant-stress exper- 
iments (creep tests). The constitutive parameters are a = 1, ai = 1, and D = 0.4. The initial 
(dimensionless) density of STZ's is Aj„, = 0.1; and the sequence of applied stresses do is as shown. 
The mesh size for numerical discretization is = 0.05. Stresses ctq are normalized so that the 
ultimate yield stress is ay = 1.132. Strains are in units eo, as defined in Eg. ( p. 12 ). Times t are in 
units Rq^, as defined in Eq.(|2.3D. 
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FIG. 4. Strain rate q, as a function of position ^, corresponding to the point t = 100 on the 
curve for ctq = 1.0 in Fig. 3. q is in units Eq Rq. 
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FIG. 5. Stress-strain curves for a sequence of simulated constant strain-rate experiments. All 
material parameters are the same as in Fig. 3. The strain rates qo are as shown. 



22 



0.02 



0.015 



q 0.01 



0.005 



3 4 



0.2 



0.15 



0.1 



0.05 



1 2 3 4 5 6 7 

FIG. 6. Two graphs of the strain rate g as a function of position ^ showing the emergence 
of shear bands near a weak spot in an otherwise stable system. The constitutive parameters are 
a = 10, o"! = 1, and P = 0.4. The uniform apphed stress is (Tq = 0.5. Initiahy, Aj^ = 0.5, A = 
everywhere except in a very narrow region near ^ = 3.5 where, at t = 0, Aj„ = 2.0, A = — 1. In (a), 
where t = 200, the system as a whole is still deforming in response to the applied stress, and the 
central shear bands arc just beginning to appear. In (b), where t = 5000, all of the deformation 
rate is concentrated in the bands, and the whole system has reached its steady state configuration. 
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